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ABSTRACT 

We present an extensive study of the double (3 model for the X-ray surface brightness 
profiles of clusters, and derive analytically the gas density and total masses of clusters 
under the hydrostatic equilibrium hypothesis. It is shown that the employment of the 
double (3 model instead of the conventional single [3 model can significantly improve 
the goodness of fit to the observed X-ray surface brightness profiles of clusters, which 
will in turn lead to a better determination of the gas and total mass distributions 
in clusters. In particular, the observationally fitted (3 parameter for the extended 
component in a double (3 model may become larger. This opens a new possibility of 
resolving the long-standing (3 discrepancy for clusters. Using an ensemble of 33 ROSAT 
PSPC observed clusters drawn from the Mohr, Mathiesen & Evrard (1999) sample, 
we find that the asymptotic value of (3fu is 0.83 ± 0.33 at large radii, consistent with 
both the average spectroscopic parameter (3 spec — 0.78 ± 0.37 and the result given by 
numerical simulations. 
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1 INTRODUCTION 

Since the pioneering work of Cavaliere & Fusco-Femiano 
(1976), the (3 model has been widely adopted in the fitting 
of the X-ray surface brightness profiles of galaxy clusters. 
In particular, despite the difference in the core radius there 
is a striking similarity between the X-ray surface brightness 
distribution predicted by the universal density profile as the 
underlying gravitational potential of a cluster and the con- 
ventional (3 model (Makino, Sasaki & Suto 1998). Yet, it has 
been realized for many years that a single (3 model is inac- 
curate and also not self-consistent. First, it fails to repre- 
sent the central excess emission associated with the cooling 
flows seen in many clusters. Second, the X-ray luminosity 
is divergent for (3 < 0.5 so that an arbitrary cutoff radius 
should be introduced (e.g. Henry & Henriksen 1986). Third, 
an increasing temperature with cluster radius is required for 
(3 < 2/3 if the asymptotic baryon fraction of cluster at large 
radii should reach a universal value defined by the Big Bang 
Nucleosynthesis (Wu & Xue 2000). Finally, the well-known 
P discrepancy (e.g. Bahcall & Lubin 1994) may also arise 
from the employment of a single (3 model. Namely, the slope 
of gas radial profile described by a single j3 model is sys- 
tematically smaller than that required by the equipartition 
between specific kinetic energy in galaxies and that in gas. 

Regardless of its simplification, the second (3 model has 
been formally adopted by a number of authors in recent 
years to represent the excess X-ray emission in the central 
cores of cooling flow clusters (e.g. Ikebe et al. 1996; 1999; 



Xu et al. 1998; Mohr, Mathiesen & Evrard 1999, MME here- 
after). An immediate consequence of fitting the X-ray sur- 
face brightness profile with a double (3 model instead of a 
single (3 model is that the resultant (3 parameters in the 
two components may both become larger. In particular, it 
seems that the (3 parameter for the extended component 
can often be greater than 2/3. Therefore, there is a possibil- 
ity that a double (3 model may allow one to resolve all the 
above puzzles associated with a single (3 model. Moreover, 
in the conventional single /3 model fit, some of the observed 
data points in the central regions of the cooling flow clus- 
ters should be omitted in order to obtain an acceptable fit 
(e.g. Jones & Forman 1984), while in some cases, the sizes 
of the excluded regions are not at all obvious. The introduc- 
tion of a double (3 model thus provides a way to quantita- 
tively describe the central excess emission. In fact, as will 
be shown in this paper, our understanding of some prop- 
erties of clusters such as the determinations of gas density 
and dynamical mass is closely connected to the issue as to 
whether or not the central narrow component is properly 
taken into account. On the other hand, there have been in- 
creasing studies both theoretically and observationally on 
the possibility that intracluster gas is multiphase (Nulsen 
1986; White & Fabian 1995; Gunn & Thomas 1996; Nagai, 
Sulkanen & Evrard 2000; Buote 2000, etc.), among which the 
detection of excess of low-energy photons in X-ray spectra 
relative to the average X-ray temperature in some clusters 
gives a convincing support to the presence of a cold gas com- 
ponent either concentrated in the central core or distributed 
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over the entire cluster. A double (3 model may correspond 
to the simplest case for multiphase medium, and therefore a 
detailed study of its properties will be helpful for our further 
investigation of the complex and multiphase intracluster gas. 

In this paper, we present the gas distribution inferred 
from a double f3 model (section 2.1) and the total cluster 
mass under the assumption that the two-phase gas is in 
hydrostatic equilibrium with the underlying gravitational 
potential of the cluster (section 2.2). We demonstrate the 
difference in the cluster properties characterized by a single 
and double (5 model using two strong cooling flow clusters, 
A2597 and A2390 (section 3). Because the sharp peak in 
the central X-ray emission can be quantitatively described 
by the second (3 model, we will be able to estimate more 
accurately the cluster mass enclosed within the central core 
and compare with that derived from strong gravitational 
lensing. We will examine whether the double /3 model cor- 
rected X-ray mass can be reconciled with the strong lensing 
result (e.g. Wu & Fang 1997). Finally, as a consequence of 
adopting a double /3 model fit, the (3 parameters in the two 
components may both become larger than the value in a 
single f3 model fit. We will study the possibility of resolving 
the well-known /3 discrepancy, using a sample of 33 clusters 
drawn from the recent work of MME. Throughout this paper 
we assume Ho = 50 km s _1 Mpc -1 and fio = 1- 



2 THE MODEL 

2.1 GAS DISTRIBUTION 

If we assume that the intracluster gas has two-phases with 
different electron temperatures Ti, where and also hereafter 
the subscript number (i = 1, 2) refers to the two-phases, the 
X-ray emission per unit volume in an energy band from E\ 
to E2, according to thermal bremsstrahlung, is 

dL x 
~dV 



(1) 



where 



3m e hc 2 \3m e c 2 / 



and 



fE 2 
JE 1 



■'/'= / 9ff(J^)e-^d^— r 



(2) 



(3) 



in which n e is the electron number density, Nz, z is the ion 
number density with atomic number Z and effective charges 
z, and is the Guant factor of the free-free emission. We 
define the average electron weight fi ei (or fj, e ) as 



(leiBei = flede = ^ N Z ,zZ 

where n e = n ei is the total electron number density, so 



(4) 



that 

E- = - 



(5) 



In particular, if intracluster gas is assumed to have the cos- 
mic primordial abundances of hydrogen and helium, then 



/j, e — 2/(1 + X) with X being the primordial hydrogen mass 
fraction X = 0.768. 

Suppose that the two-phase emission gives rise to, cor- 
respondingly, the two components in the observed X-ray sur- 
face brightness profile described by a double /3 model: 

/ 2 \ -3/3,+O.B 

S(r) = 5>«( 1 + ^-) . (6) 



By inverting S(r) and comparing with eq.(l), we have (e.g 
Cowie, Henriksen & Mushotzky 1987) 

4^ 2 



(OKi(o) 

and 
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Heinli = Hei(0)nli(0) ( 1 + — 



-30i 



(7) 



(8) 



Eliminating fi e i in the above equations in terms of eqs.(4) 
and (5), we can get the electron number densities for the 
two components as well as the combined electron number 
density n e 



n e i(0) 



1 + 



-I 1/2 



-3ft 



-3ft 



,(9) 



or 

_ f n eO\ ~ 

\ n e J 

and 



in which 



h ei = Tlei(O) 1 + 



1/2 



-3ft 



(10) 



(11) 



(12) 



Note that although eq.(ll) can be formally written as n gas — 
[ngasi + r *gas2] 1//2 ) which has been adopted in some recent 
work (e.g. Reiprich & Bohringer 1999), the determination of 
the central gas density for each component (n 9as i or n ga3 2) 
described by a single j3 model is dependent on both compo- 
nents of S(r). In terms of eq.(7), the central electron number 
density, n e ;(0), is related to the observationally fitted central 
surface brightness Soi through 



4tt 1/2 r(3ft) S 0l 
a{T t ) gil i e T(3/3 l - |) r ci 



nlm = -^r, (13) 



in which 

1 _ 1 s 0j r ci g s (Ti\ 1/2 TO r(3A - i) 

An r,:iHj \T :i J r(3^) T(3/3 i ' 



'1 ~ D' 



where and also hereafter j — 1, 2 and j 7^ i. Of course, one 
can also estimate n e j(0) using the total X-ray luminosity L x 
of the cluster (eq.[l]). For example, in the case of isothermal 
gas distributions for the two-phase gas we have 



nlm~ L ' A " B <- 



r(3A) 



7 r3/2a(T i ) 5< Mer- c 3 i T(3ft - |) : 



(15) 
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where 



B 



■ 1 Scy Ki 3ft 



So; r c 2 ; 3ft - | 



(16) 



Apparently, the total X-ray luminosity can be estimated 
simply from 



E 



4-iT 2 r 2 ci Soi 

3 
2 



3ft 



(17) 



2.2 CLUSTER MASS 

The total mass in gas within a sphere of radius r is 

1 1/2 



M 9as (r) 



: 47r^ e m„n e0 / — — 

io V 



n e i(r) 



r 2 dr, 



(18) 



If the X-ray emitting gas is in hydrostatic equilibrium with 
the underlying gravitational potential of the cluster, the to- 
tal dynamical mass within r is 

„2 



M(r) 



v-* d(riikTi) 



Gfj, e m p n e ^—i dr 



(19) 



in which n» is the total particle number density for the ith- 
phase gas. If the abundances of hydrogen and helium for a 
single-phase gas are assumed to remain unchanged over the 
whole cluster, we have where p s ; is a propor- 

tionality coefficient. If the two-phases both have the cosmic 
mixed-abundances of hydrogen and helium, // s j = 1.934. 
Next, we assume a poly tropic equation of state for each 
single-phase gas, namely, T, = Tio[ne;/n e ;(0)] 7i_1 . In this 
case, eq.(19) can be written to be 



M{r) 



E 



7ifcT i0 r / n 



GiMrrip \ n ei (0) 



where [n — fie/fJ-si- A straightforward computation using 
the electron number density found in the above subsection 
yields 



M(r)= ^M.(r)g) 



1 + 11 i r +<i 



(21) 



where M»(r) is the total cluster mass determined by assum- 
ing that the X-ray emission is from a single-phase medium 
described by a single /3 model: 

-i 



M t {r) = 3ft 7i 



n e i(Q) 



(22) 



It is easy to show that eq.(21) reduces to the form of eq.(22) 
if the two-phase gas components are assumed to be identical 
or the second component vanishes. 



3 APPLICATION TO X-RAY CLUSTERS 
3.1 A2579 

We first choose a typical cooling flow cluster, A2579, to 
demonstrate how a double (3 model works. A2579 has been 
extensively studied by Sarazin & McNamara (1997). They 



showed that the merged ROSAT HRI and PSPC surface 
brightness profile S(x) cannot be fitted by a single /3 model 
unless the data points within the central region of radius 
of 0.18 Mpc are removed. Meanwhile, the core radius has 
not been well constrained due to the presence of the cen- 
trally peaked X-ray emission, and the gas and total masses 
within the central region derived from the deconvolved gas 
density are apparently larger than those obtained from the 
best-fit single /3 model. We now make an attempt to fit 
the same data set S(x) using a double j3 model (eq.[6]). 
Fig.l displays the observed and our best- fit X-ray surface 
brightness profiles of A2579. It turns out that the good- 
ness of the fit has been significantly improved when a dou- 
ble (5 model is applied to the entire data points, yielding 
X 2 /v = 41.16/26. Recall that the minimum x 2 /v is 59.99 
for 29 d.o.f in a single (3 model fitting (Sarazin & McNa- 
mara 1997). The best-fit values of the (3 parameter and 
core radius for the narrow and extended components are 
(ft,r c i) = (0.70 ± 0.32, 0.047 ±0.016 Mpc) and (ft>,r c2 ) = 
(0.66±0.03, 0.15±0.05 Mpc), respectively. Note that our fit- 
ted f3 parameter for the extended component is close to the 
value (0.64) obtained by Sarazin & McNamara (1997) who 
have excluded the central cooling flow region. Alternatively, 
the two components in our double j3 model exhibit similar f3 
values. The same situation has been known for A1795 (Xu 
et al. 1998), while MME assumed the two components of 
a double /? model to have the same j3 in their analysis of 
45 nearby clusters. We have also tried the double /3 fit by 
correcting the PSF of the ROSAT HRI with a FWHM of 4 
arcseconds. This results in a slightly smaller core radius for 
the narrow component (ft,r c i) = (0.67 ± 0.33, 0.041 ± 0.016 
Mpc), while the extended component remains almost un- 
changed. Finally, we have compared in Fig.l the electron 
number density given by the deprojection technique with 
that calculated from our best-fit double f3 model assuming 
Ti = 2.0 kcV and T2 = 4.34 keV for the narrow and extended 
components, respectively. The agreement between the two 
results is good with \ 2 = 41.06 for 32 d.o.f. although the 
very central point shows a deviation from the expectation 
of our double j3 model. Note that the temperature of the 
cluster, especially the outer cluster, is poorly constrained. 
Our choice of 7\ and T2 is simply based on the single tem- 
perature model fits by Sarazin & McNamara (1997) to the 
ROSAT PSPC spectra of the core region (0.1 Mpc) and the 
outer cluster (0.25-1 Mpc), respectively. 

3.2 A2390 

The reason why we choose A2390 for a further demonstra- 
tion of the double (5 model is as follows: First, a single /3 
model fails to fit the merged ROSAT HRI and PSPC sur- 
face brightness profile S(x) of the cluster, which has been 
well observed out to a radius of ~ 2 Mpc (Bohringer et 
al. 1998); Second, the mass distribution of the cluster has 
been mapped with the weak lensing technique (Squires et 
al. 1996). It deserves to be examined whether the cluster 
mass derived from a double /3 model can be reconciled with 
the weak lensing result. The latter is believed to provide a 
reliable mass estimate, independently of the cluster matter 
contents and their dynamical states. Third, none of the sug- 
gested models in literature based on the fitting of S x have 
reproduced the mass profile given by the weak lensing analy- 
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sis (Bohringer et al. 1998). In Fig. 2 we plot the X-ray surface 
brightness profile of A2390 observed from the ROSAT HRI 
and PSPC (Bohringer et al. 1998), along with a single and 
double f3 model fitting. The best-fit parameters are (/?, r c ) = 
(0.56 ± 0.01,0.11 ± 0.01 Mpc) for a single [3 model with 
X 2 /y = 209.08/40, and (/3i,r c i) = (0.92 ± 0.22, 0.46 ± 0.16 
Mpc) and ((3 2 ,r c2 ) = (0.51 ± 0.01,0.046 ± 0.010 Mpc) for 
a double (3 model with x 2 l v = 52.48/37, respectively. In- 
deed, a single j3 model fit is not acceptable. Bohringer et 
al. (1998) obtained a good fit by introducing a Gaussian 
peak plus a single f3 model. Our double (3 model does 
provide a significantly reduced \ 2 fit to the entire data 
points although the two components appear to be unusual 
in shape. Because A2390 is located at moderate redshift 
(z = 0.228), we may need to consider whether the core ra- 
dius is overestimated in the above fitting due to the PSF 
of the ROSAT HRI. We repeat the double /3 model fitting 
using a PSF with a FWHM of 4 arcseconds. The best fit pa- 
rameters become (f3i,r cl ) = (0.98 ± 0.12, 0.49 ± 0.08 Mpc) 
and C32,r c2 ) = (0.50 ±0.01, 0.035 ±0.005 Mpc). Indeed, for 
the second component the core radius is apparently reduced, 
while the (3 parameter remains roughly the same. 

Now, we calculate the projected cluster mass within ra- 
dius r. The major uncertainty here is the X-ray tempera- 
tures for the two components of the double (3 model. The 
presence of the cooling flow may lead to an underestimate of 
the overall temperature of the cluster. The reported temper- 
ature of the cluster in the literature ranges from the average 
value 7.7lJ;gg keV over the cluster (Rizza et al. 1998) to 
the cooling flow corrected value 23.09 ± 35.12 keV (White 
2000). Moreover, there is an apparent discrepancy between 
the ROSAT PSPC and ASCA CIS measured temperatures 
within the central region of the cluster. In the following we 
first adopt the two-temperature model of Bohringer et al. 
(1998), which consists of a narrow, low-temperature compo- 
nent with Ti and an extended, high-temperature one with 
T2. The projected cluster mass deriving from this double 
(3 model based on the best-fit values of Ti = 2.0 keV and 
T 2 = 11.5 keV by Bohringer et al. (1998) is shown in Fig.2. 
While the shape of the resulting mass profile seems to be 
consistent with the weak lensing result, the double f3 model 
with the above temperature data has significantly underes- 
timated the cluster mass as compared with the weak lensing 
analysis. Next, we treat the two temperatures Ti and T2 
as free parameters and search for (Ti, T2) that give rise to 
the best fit of the derived cluster mass from the double (3 
model to the weak lensing result. This yields Ti = 4.32 keV 
and T2 = 20.0 keV for the narrow and extended compo- 
nents, respectively, with a reduced \ 2 jv value of 2.38. Inter- 
estingly, these best-fit temperatures for the narrow and ex- 
tended components are consistent with the reported value of 
~ 4 keV in the central region of the cluster inside ~ 0.5 Mpc 
(Bohringer et al. 1998) and the cooling flow corrected value 
23.09 ± 35.12 keV (White 2000), respectively. We have also 
displayed in Fig.2 the projected mass profiles derived from 
other two models for S x '- (l)a single (3 model and (2)a single 
(3 model plus a Gaussian peak. For these two models we use 
an isothermal plasma of T — 9.0 keV, a value obtained from 
the ASCA spectra within a radius of 1.7 Mpc (Bohringer et 
al. 1998). The goodness of the fits of the resulting cluster 
masses to the weak lensing measurement reads X 2 l v = 3.55 
and 4.63 for models (1) and (2), respectively. In a similar 



way to the above analysis for the double f3 model, we have 
performed the x fits of the projected cluster masses given 
by models (1) and (2) to the weak lensing determined mass 
profile by treating the temperature as a free parameter. It 
turns out that the best-fit values of temperature for the two 
models are T = 8.74 keV with xV" = 3 - 48 and T = 9 - 15 
keV with x 2 l v = 4.59, respectively. These values are well 
within the ROSAT and ASCA spectral measurements. It 
is likely that the double (3 model fit gives the smallest x 2 
value, and thus provides the most precise description for 
the intracluster gas, although the current fit is still unsat- 
isfactory. Nevertheless, it should be pointed out that the 
measurement uncertainties in the fitting of S x and the spec- 
tral measurement of temperature have not been included in 
the final mass estimate of the cluster. In particular, using a 
double (3 model instead of a single (3 model will lead to a 
large uncertainty in the derived cluster mass since the cen- 
tral surface brightness, core radius and f3 parameter in a 
double (3 model fit cannot be well constrained at present. 
Finally, as compared with the single (3 model, the double 
(3 model results in a significantly large cluster mass within 
small radii (~ 100 kpc), and the overall dark matter dis- 
tribution clearly exhibits two distinct length scales. Similar 
result has also been reported in previous studies (e.g. Ikebe 
et al. 1996). However, it is unlikely that the central mass 
excess as a result of the employment of the double (3 model 
can resolve the mass discrepancy between the strong lens- 
ing and X-ray measurements (Wu & Fang 1997), despite the 
fact that our fitted core radius (r C 2 = 0.035 ± 0.005 Mpc) is 
relatively small. Our result agrees with the recent analysis 
by Lewis et al. (1999), who found the X-ray mass within the 
arc radius is 2.3 times smaller than the value reported by 
Allen (1998), 2.1 x 10 14 M Q . This indicates again that the 
previously claimed mass discrepancy between strong lensing 
and other methods has probably arisen from the oversim- 
plification of lensing model for matter distribution in the 
central cores of the strong lensing clusters. 



4 THE 13 DISCREPANCY 

In this section we discuss whether the well-known puzzle, 
the so-called (3 discrepancy (e.g. Bahcall & Lubin 1994), 
may have also arisen from the employment of the conven- 
tional (3 model for intracluster gas. Briefly, if both galaxies 
and gas are the tracers of the depth and shape of a common 
gravitational potential of a cluster, we would expect the fol- 
lowing identity (Cavaliere & Fusco-Femiano 1976; Bahcall 
& Lubin 1994) 



f3spec — 



a 2 d\nngas/d\nr 
kT/fmip d\nn ga i/d\nr 



= Pfit, 



(23) 



where a is the velocity dispersion of galaxies, and the 
number densities of galaxies and gas are denoted by n ga i 
and n gas , respectively. For simplicity, we have assumed an 
isotropic and isothermal matter distribution for the clus- 
ter. The (3 discrepancy arises from a statistical estimate of 
f3 spec and f3f it , which turns out that on average f3 Bp ec (~ 1) 
is larger than f3ju (~ 2/3). Basically, two solutions to the 
puzzle have so far been suggested: Either the galaxy pro- 
file in clusters may be less steeper than the simplified King 
model, or the present X-ray and optical observations have 
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only probed some finite regions of clusters so that the asymp- 
totic value of (3 fu has not yet been reached (Gerbal, Durret 
& Lachieze-Rey 1994). 

Numerical studies of cluster formation and evolution 
(e.g. Navarro, Frenk & White 1995; Lewis et al. 1999) 
have shown that the gas density profiles of clusters become 
steeper with radius and therefore, a larger (3 value can be ob- 
tained if the fit to the observed X-ray surface brightness pro- 
file is extended to large radii. Indeed, Vikhlinin et al.(1999) 
have tried a single f3 model fit to the gas distribution in 
the outer regions of clusters between 0.3 and 1 virial radius. 
They found that the typical (3 values among the 39 clusters 
observed by ROSAT PSPC would be larger by « 0.05 than 
that derived from the global fit. These results provide a use- 
ful clue to our reconsideration of the f3 discrepancy because 
only the asymptotic /3 value at large radii is of significance 
for the evaluation of (3 fa- As has been pointed out in the 
above section, the f3 parameters for the two components of 
a double (3 model may both become larger than the value 
in a single (3 model fit. Replacing the (3 in eq.(23) by the (3 
parameter of the extended component of a double (3 model 
may give rise to a larger value of f3fu- Hence, this opens a 
new possibility of reconciling f3 apec with (3fu- 

4.1 CLUSTER SAMPLE 

We begin with the sample of 45 ROSAT PSPC observed 
nearby clusters compiled recently by MME. We choose this 
sample because MME have already tried a double (3 model 
fit to the 18 clusters of their sample, which can be used for 
the purpose of comparison, although their double j3 model 
assumed the same /3 parameter for the two components. In 
order to estimate (3 ape c , we exclude those clusters whose ve- 
locity dispersion and temperature are not observationally 
determined. This reduces the MME sample to 33 clusters 
(see Table 1 and Table 2). We reanalyze the archival ROSAT 
PSPC imaging data of the 33 clusters following the same 
procedure of MME, who binned the X-ray surface brightness 
according to photon counts. We have also tried to measure 
S x (r) using the concentric rings of equal width (1 arcminute) 
but found that the data points are too few for some clusters. 
To facilitate the double (3 model fit to S x (r), we will use the 
results of S x (r) from MME. 

We fit the observed X-ray surface brightness of clusters 
to a double f3 model described by eq.(6). We employ the 
Monte Carlo simulations and the ^-Ht to obtain the best- 
fit parameters (So;, r ci , f3i) and estimate their error bars, 
in which we keep the same outer radii of the fitting regions 
as those defined by MME. Essentially, the fitting can be 
classified as two types, and a typical example for each type 
is displayed in Fig. 3. In the first case which is applied to 
15 clusters in Table 1, the narrow component contributes 
mainly to the central X-ray emission and exhibits a sharp 
drop at large radii, while the extended component domi- 
nates the emission in the outer region. As for our example 
A2052 in Fig.3, a single (3 model fit gives (3 = 0.593 ± 0.050 
but the minimum \ 2 IS 310.67 for 49 d.o.f. Our double (3 
fit yields f3i = 0.810 ± 0.194 and (3 2 = 0.773 ± 0.029 with 
X 2 jv = 60.67/46. It appears that our best-fit (3 parameter 
for the extended component is comparable to the one (0.712) 
found by MME using a double [3 model fit with the same (3 
parameter for the two components. In the second case, the 



central emission is similarly governed by the narrow com- 
ponent. However, the narrow component exhibits a flatter 
slope than the extended component does in the outer region. 
Yet, the amplitude of the narrow component is significantly 
smaller than that of the extended component. We will thus 
take a less vigorous approach to the selection of the (3 pa- 
rameter: The value of the extended component will be used 
in our evaluation of (3f it below. Mathematically, the fitted 
X-ray surface brightness profile will be eventually dominated 
by the narrow and shallower component at very large radii 
because of its smaller value of (3. For A1060 shown in Fig.3, 
the minimum \ 2 °f 598.3 in a single (3 model fit for 113 
d.o.f has been reduced to 139.9 in our double [3 model fit for 
110 d.o.f. The best-fit [3 parameter of the extended compo- 
nent is 0.732 ± 0.065, in comparison with the MME result 

;9 = 0.7031°;°**- 

The best-fit f3 parameters by a double (3 model for our 
cluster sample are listed in Table 1, together with the results 
for a single (3 model fit, Quoted errors are 68% confidence 
limits. Note that the (3 values of the narrow components for 
about 1/3 of our clusters have not been well constrained. 
This is partially due to the fact that a single (3 model fit is ac- 
ceptable for some clusters such as A401, A1367, A1651 and 
A4059. Another reason comes from the sparse data points 
at small radii, which makes it difficult to place a robust con- 
straint on the narrow component. 

We searched the literature for the velocity dispersion of 
galaxies and the temperature of X-ray emitting gas in our 
cluster sample (Table 2). Essentially, we take the X-ray tem- 
perature data of White (2000) (except A2244 and A3532) 
because the effect of the cooling flow on the cluster temper- 
ature has been corrected for. The majority of the velocity 
dispersion data are chosen from the substructure corrected 
values by Girardi et al. (1998). The f3 spec value (the signif- 
icance at 68% confidence) is correspondingly calculated for 
each cluster and listed in Table 2. 

4.2 THE (3 PARAMETERS 

The mean value of [3 spec over 33 clusters in Table 2 is 
(flspec) = 0.78 ±0.37, which is consistent with previous find- 
ings within uncertainty [see Wu, Fang & Xu (1998) for a 
summary]. This value can be slightly raised to (f3 spec ) — 
0.80 ± 0.36 if A754 is excluded from the list. This cluster 
gives rise to the smallest value of j3 apec because of its rela- 
tively small velocity dispersion or high temperature. Simi- 
lar effect was also noticed in previous work, i.e., the mean 
value of f3 sp ec can be biased high or low because of a few 
clusters in the sample. For instance, Bird, Mushotzky & 
Metzler (1995) found that their mean value of (3 sp ec = 1-20 
drops to 1.09 if A2052 is removed from their cluster list. Us- 
ing a large sample of 149 clusters, Wu, Fang & Xu (1998) 
showed (/3 spec ) = 1.00 ± 0.49 but the median value is only 
(f3 spec ) = 0.80l ;°*. 

The mean values of f3fu (we will assume a simplified 
King model for galaxy number density profile) from Table 1 
are {(3 fit ) = 0.61 ± 0.09 and ((3 fit ) = 0.83 ± 0.33 for a single 
and double f3 fit, respectively. While the error bar is still 
large, the value of (/3/;t) has been apparently increased in 
the case of a double [3 model fit to S x (r). For comparison, the 
mean value of (3f it in the MME sample is {f3fu) = 0.61 for 
a single f3 model fit among 27 clusters, while (/?/;«) = 0.70 
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Table 1. Cluster Sample: f3fu 
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for a double /3 model fit for the rest 18 clusters. Recall that 
their double f3 model assumed the same fi parameter for the 
two components. Our result of f3f it for a single [3 model is in 
good agreement with the one found by MME. However, our 
double f3 model fit yields a slightly larger mean value of /9/it, 
which is nevertheless consistent with the claim (0.84 ± 0.1) 
by Bahcall & Lubin (1994) (Note the different definition of 
f3fu), and in particular with the numerical result (0.82 ± 
0.06) by Navarro et al. (1995). 

We now examine whether the observationally fitted pa- 
rameter (3 f it in the outer regions of clusters have any depen- 
dence on temperature or velocity dispersion. Fig. 4 shows 
the best-fit f3f it versus T for our cluster sample. We fit 
the data points to an expression f3fu = (T/keV) a + b, 
and find (a,b) = (0.17 ± 0.03, -0.73 ± 0.05) and (a, b) = 
(0.22 ± 0.03, -0.71 ± 0.08) for a single and double [3 model 
fit, respectively. We thus confirm the claim by Vikhlinin et 



al. (1999) that (3 fa has a mild trend to increase with tem- 
perature, although we get somewhat a larger variation of 
(3 fi t from f3 fit « 0.56 (0.48) for T = 3 keV to w 0.93 (0.75) 
for T = 10 keV in a double (single) model fit. 

We then study the distribution of f3 spec /f3f it for each 
cluster, aiming at examining whether the spectroscopic pa- 
rameter f3 spec correlates with the observationally fitted pa- 
rameter (3 fu from S x (r). To this end, we display in Fig.5 
[3 S pec against (3fu for our 33 clusters. Indeed, there is 
a positive correlation between the two (3 parameters. A 
X 2 fit excluding the data point of A754 yields [3 spec — 

10 p.018±0.044 / gl.ia±0.29 ) indicating that a high corre _ 

sponds to a high (3 3pec . Moreover, the above relation is also 
consistent with f3 spec — (3f it , although the mean values of 
f3 spec and (3f it are on average smaller than unity (see be- 
low). 

Finally, we examine whether there is a (3 discrepancy 
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Table 2. Cluster Sample: f3 sp ec 
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References.- (a)Girardi et al. 1998; (b)Whitc 2000; (c)Zabludoff et 
al. 1990; (d)Jones & Forman 1999; (e)Bird et al. 1995; (f)Struble 
& Rood 1991; (g)Mushotzky & Scharf 1997; (h)Fadda et al. 1996; 
(i)Edge et al. 1990; (j)Green et al. 1988; (k)Owen et al. 1997. 



for clusters, we plot in Fig. 6 the ratio of f3 spec to f3f it for 
each cluster in our cluster sample. For a single j3 model, the 
mean ratio is (f3 spec / /3fit) = 1-28 ± 0.41, while for a dou- 
ble (3 model, this value becomes {f3 sp ec/ ' Pfit) = 0.99 ± 0.35. 
Therefore, it seems that there is no apparent (3 discrepancy 
for clusters within the framework of a double /3 model fit to 
S x (r). 



5 DISCUSSION AND CONCLUSIONS 

The conventional (3 model appears to be inaccurate, and 
the extrapolation of a single /3 model to large radii could be 
misleading. The monotonic increase of baryon fraction with 
radius and the so-called /3 discrepancy, together with the 



luminosity divergence as a result of the too small (3fu pa- 
rameters {f3 fu < 1/2) for some clusters, may have a common 
origin: A single (3 model fit to the entire regions of the X- 
ray surface brightness profiles of clusters yields an underes- 
timate of j3 fu parameter at large radii. The exclusion of the 
central data points or the employment of a double f3 model 
in the fitting of the observed X-ray surface brightness pro- 
files can lead to an apparent increase of j3fu at large radii, 
which may allow one to resolve all the puzzles mentioned 
above, although these two empirical methods both depend 
on the extent of the fitting regions. Yet, in the conventional 
treatment the cooling flow regions are often excluded from 
the (3 model fit, which may suffer from the uncertainty of 
different strategies of choice of the cutoff radii (e.g. Jones 
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k, Forman 1984; White, Jones & Forman 1997; Vikhlinin, 
Forman & Jones 1999). As a result, the accuracy of our 
determination of the total cluster mass from X-ray mea- 
surements may be correspondingly affected. Employing the 
second /3 model to quantitatively describe the central excess 
X-ray emission can significantly improve the goodness of the 
fit, which thus enables us to derive more accurately the gas 
and total dynamical masses of clusters. In the present pa- 
per, we have extensively studied the properties of a double 
f3 model, and demonstrated its applications using two ex- 
amples: A2597 and A2390. In particular, as compared with 
a single /3 model, a double (5 model via the hydrostatic equi- 
librium hypothesis for intracluster gas gives rise to a larger 
X-ray cluster mass at the central region, although this cor- 
rection is still insufficient to resolve the discrepancy between 
the strong lensing and X-ray determined cluster masses en- 
closed within the arc radius (e.g. Wu et al. 1998). It is hoped 
that the mathematical treatments of a double /3 model de- 
veloped in the present paper will be useful for future work. 

The (3 parameter as an indicator of the dynamical prop- 
erties of clusters discussed in the present paper and the 
baryon fraction of clusters as an indicator of cosmologi- 
cal matter composition discussed in literature (White et al. 
1993; David, Jones & Forman 1995; Wu & Xue 2000 and 
references therein) concern mainly the asymptotic behavior 
of clusters at large radii, while a double (3 model provides a 
more precise way to describe the variation of intracluster gas 
at large radii. We have thus studied the possibility of resolv- 
ing the well-known /3 discrepancy using a double model fit 
to the X-ray surface brightness of clusters. We have found, 
based a sample of 33 clusters drawn from a recent catalog 
of MME, that the spectroscopic parameter /3 apec is on aver- 
age consistent with the observationally fitted parameter /3/it 
from the X-ray surface brightness profile. 

Of course, the definite resolution to the /3 discrepancy 
also needs a better understanding of the distributions of 
cluster galaxies and their velocity dispersion at large radii 
(Bahcall & Lubin 1994; Girardi et al. 1996), for which we 
have assumed a simplified King model with a constant and 
isotropic velocity dispersion in the present paper. A shal- 
lower galaxy density profile (e.g. n ga i tx r~ 2A ) can give 
rise to a higher value of /3/it- Alternatively, the second X- 
ray core radius will become larger in a double f3 model fit. 
In our cluster sample, the mean core radii of the narrow 
and extended components are (r c i) = 0.14 ± 0.14 Mpc and 
(r c i) = 0.48 ± 0.38 Mpc, respectively. Since the extended 
component is involved in the evaluation of /3f it , a good ap- 
proximation of /3/ it can only be achieved at very large radii. 
So, the effect of finite observational sizes may cause another 
uncertainty, and the extension of the fitting region will lead 
to an increase of /3/it- As a whole, these two effects, the 
shallower galaxy density profile and the extension of the 
X-ray fitting region, can both result in a high /3f it . There- 
fore, it is possible that the /3 discrepancy still remains but 

(3 fit ^* /3spec- 

Finally, within the framework of energy equipartition 
the ratio of kinetic energy in the galaxies to that in the 
gas should be close to unity. Although our estimate of 
{Pspec) = 0.78±0.37 is still consistent with such a prediction 
within the error bars, this value, together with the consis- 
tent /3/it parameter (/3/it) = 0.83 ± 0.33, does not exclude 
the possibility that the actual value of (/3 spe c) is smaller than 



1, namely, (a 2 ) < {kT} / [im p . Here, we will not intend to 
explore in detail the physical mechanisms for the scenario. 
However, we would like to point out that the dynamical fric- 
tion would resist the motion of galaxies in clusters, leading 
to a velocity dispersion smaller than that expected from the 
purely dissipationless hypothesis (e.g. Carlberg & Dubinski 
1991; Metzler & Evrard 1994). This may account for our 
result - why galaxies appear to be cooler than intracluster 
gas. 
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r (Mpc) 



Figure 1. Upper panel: The merged ROSAT HRI and PSPC surface brightness profile for A2579 (Sarazin <fc McNamara 1997), along 
with the best-fit double /3 model (solid line). The effect of the PSF of the ROSAT HRI has been corrected. Lower panel: A comparison 
of the electron number density derived from the deprojection technique (Sarazin & McNamara 1997) and that from the best-fit double 
/3 model (solid line). The dotted lines represent the two components of the double f3 model. 
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Figure 2. Upper panel: The superimposed ROSAT HRI (squares) and PSPC (circles) surface brightness profiles for A2390 (Bohringcr 
et al. 1998). Also plotted are the best-fit double fi model (solid line) and the best-fit single /3 model (dashed line) for the entire data 
points. The dotted lines represent the two components of the double (i model. Lower panel: A comparison of the projected cluster masses 
revealed by weak gravitational lensing (filled circles) (Squires et al. 1996) and by the X-ray measurements via hydrostatic equilibrium 
hypothesis: thin solid line - the best-fit double (i model with Ti = 2 keV and T2 = 11.5 keV; thick solid line — the best-fit double (3 model 
with Ti = 4.32 keV and T2 = 20.0 keV; dashed line - the best-fit single (3 model for the entire data points; dotted line - the single /3 
model fit by exclusion of the central emission excess. The strong lensing derived cluster mass is also shown by asterisk (Wu et al. 1998). 



© 0000 RAS, MNRAS 000, 000-000 



12 Xue and Wu 




Figure 3. Two typical examples of the ROSAT PSPC observed and the double fi model fitted surface brightness profiles, A2052 and 
A1060. The dotted lines represent the two components of the double j3 model. 
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Figure 4. The /3 parameters obtained from a single (upper panel) and double (lower panel) f3 model fit plotted against the X-ray 
temperature. The dotted lines are the best-fit correlations represented by the form fi = T a + b, in which a = 0.17 ± 0.03 and 0.22 ± 0.03 
for a single and double (i model, respectively. Cygnus-A is not shown in the plot due to its high temperature (39.40 keV). 
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Figure 5. The spectroscopic parameter f5 spec is plotted against the observationally fitted parameter (3 fit from S x (r) for our sample of 
33 clusters. The solid line is the best x 2-n t power-law relation to the data set with exclusion of A754: /3 spec oc /3^ 2 , while the dotted 
line denotes f3 spec = /3fu- 
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Figure 6. The ratio of f3 spec to fiji t for an ensemble of 33 clusters. Upper and lower panels correspond to a single and double (i model 
fit, respectively. The dotted lines denote the mean ratios. 
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